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Abstract 

Signals of QCD instanton-induced processes are searched for in deep-inelastic scattering 
(DIS) at the electron-proton collider HERA in a kinematic region defined by the Bjorken- 
scaling variables x > 10~ 3 , 0.1 < y < 0.6 and photon virtualities 10 < Q 2 < 100 GeV 2 . 
Several observables characterising hadronic final state properties of QCD instanton-induced 
events are exploited to identify a potentially instanton-enriched domain. While an excess 
of events with instanton-like topology over the expectation of the standard DIS background 
is observed it can not be claimed to be significant given the uncertainty of the simulation. 
Upper limits on the cross-section for instanton-induced processes of between 60 pb and 
1000 pb are set dependent on the kinematic domain considered. The data do not exclude 
the cross-section predicted by instanton perturbation theory for small instanton sizes. At 
large instanton sizes a naive extrapolation of instanton perturbation theory yields a cross- 
section in the range of sensitivity of this study. Such a cross-section is not observed, in 
agreement with non-perturbative lattice simulations of the QCD vacuum. 



To be submitted to Eur. Phys. J. C 



C. Adloff 33 , V. Andreev 24 , B. Andrieu 27 , T. Anthonis 4 , A. Astvatsatourov 35 , A. Babaev 23 , 

J. Bahr 35 , P. Baranov 24 , E. Barrelet 28 , W. Bartel 10 , J. Becker 37 , M. Beckingham 21 , 

A. Beglarian 34 , O. Behnke 13 , C. Beier 14 , A. Belousov 24 , Ch. Berger 1 , T. Berndt 14 , J.C. Bizot 26 , 

J. Bohme 10 , V. Boudry 27 , W. Braunschweig 1 , V. Brisson 26 , H.-B. Broker 2 , D.P. Brown 10 , 

W. Bruckner 12 , D. Bruncko 16 , F.W. Biisser 11 , A. Bunyatyan 12 ' 34 , A. Burrage 18 , 

G. Buschhorn 25 , L. Bystritskaya 23 , A.J. Campbell 10 , T. Carli 10 , S. Caron 1 , 

F. Cassol-Brunner 22 , D. Clarke 5 , C. Collard 4 , J.G. Contreras 7 - 41 , Y.R. Coppens 3 , 

JA. Coughlan 5 , M.-C. Cousinou 22 , B.E. Cox 21 , G. Cozzika 9 , J. Cvach 29 , J.B. Dainton 18 , 
W.D. Dau 15 , K. Daum 33 ' 39 , M. Davidsson 20 , B. Delcourt 26 , N. Delerue 22 , R. Demirchyan 34 , 

A. De Roeck 10 - 43 , EA. De Wolf 4 , C. Diaconu 22 , J. Dingfelder 13 , P. Dixon 19 , V. Dodonov 12 , 
J.D. Dowell 3 , A. Droutskoi 23 , A. Dubak 25 , C. Duprel 2 , G. Eckerlin 10 , D. Eckstein 35 , 

V. Efremenko 23 , S. Egli 32 , R. Eichler 36 , F. Eisele 13 , E. Eisenhandler 19 , M. Ellerbrock 13 , 
E. Elsen 10 , M. Erdmann 10 ' 40 ' 6 , W. Erdmann 36 , P.J.W. Faulkner 3 , L. Favart 4 , A. Fedotov 23 , 
R. Feist 10 , J. Ferencei 10 , S. Ferron 27 , M. Fleischer 10 , P. Fleischmann 10 , Y.H. Fleming 3 , 

G. Fliigge 2 , A. Fomenko 24 , 1. Foresti 37 , J. Formanek 30 , G. Franke 10 , G. Frising 1 , 

E. Gabathuler 18 , K. Gabathuler 32 , J. Garvey 3 , J. Gassner 32 , J. Gayler 10 , R. Gerhards 10 , 

C. Gerlich 13 , S. Ghazaryan 4 - 34 , L. Goerlich 6 , N. Gogitidze 24 , C. Grab 36 , V. Grabski 34 , 

H. Grassier 2 , T. Greenshaw 18 , G. Grindhammer 25 , T. Hadig 13 , D. Haidt 10 , L. Hajduk 6 , 
J. Haller 13 , W.J. Haynes 5 , B. Heinemann 18 , G. Heinzelmann 11 , R.C.W. Henderson 17 , 
S. Hengstmann 37 , H. Henschel 35 , R. Heremans 4 , G. Herrera 7 ' 44 , 1. Herynek 29 , 

M. Hildebrandt 37 , M. Hilgers 36 , K.H. Hiller 35 , J. Hladky 29 , P. Hoting 2 , D. Hoffmann 22 , 

R. Horisberger 32 , A. Hovhannisyan 34 , S. Hurling 10 , M. Ibbotson 21 , C. Igsever 7 , M. Jacquet 26 , 

M. Jaffre 26 , L. Janauschek 25 , X. Janssen 4 , V. Jemanov 11 , L. Jonsson 20 , C. Johnson 3 , 

D. P Johnson 4 , M.A.S. Jones 18 , H. Jung 20 ' 10 , D. Kant 19 , M. Kapichine 8 , M. Karlsson 20 , 
O. Karschnick 11 , F. Keil 14 , N. Keller 37 , J. Kennedy 18 , 1.R. Kenyon 3 , S. Kermiche 22 , 

C. Kiesling 25 , P. Kjellberg 20 , M. Klein 35 , C. Kleinwort 10 , T. Kluge 1 , G. Knies 10 , B. Koblitz 25 , 
S.D. Kolya 21 , V. Korbel 10 , P. Kostka 35 , S.K. Kotelnikov 24 , R. Koutouev 12 , A. Koutov 8 , 

J. Kroseberg 37 , K. Kruger 10 , T. Kuhr 11 , T. Kurca 16 , D. Lamb 3 , M.P.J. Landon 19 , W. Lange 35 , 
T. Lastovicka 35,30 , P. Laycock 18 , E. Lebailly 26 , A. Lebedev 24 , B. LeiBner 1 , R. Lemrani 10 , 
V. Lendermann 7 , S. Levonian 10 , M. Lindstroem 20 , B. List 36 , E. Lobodzinska 10,6 , 

B. Lobodzinski 610 , A. Loginov 23 , N. Loktionova 24 , V. Lubimov 23 , S. Liiders 36 , D. Liike 710 , 
L. Lytkin 12 , N. Maiden 21 , E. Malinovski 24 , 1. Malinovski 24 , S. Mangano 36 , R. Maracek 25 , 

P. Marage 4 , J. Marks 13 , R. Marshall 21 , H.-U. Martyn 1 , J. Martyniak 6 , S.J. Maxfield 18 , 

D. Meer 36 , A. Mehta 18 , K. Meier 14 , A.B. Meyer 11 , H. Meyer 33 , J. Meyer 10 , P.-O. Meyer 2 , 
S. Mikocki 6 , D. Milstead 18 , S. Mohrdieck 11 , M.N. Mondragon 7 , F. Moreau 27 , A. Morozov 8 , 
J.V. Morris 5 , K. Miiller 37 , P. Murin 16,42 , V. Nagovizin 23 , B. Naroska 11 , J. Naumann 7 , 

Th. Naumann 35 , G. Nellen 25 , PR. Newman 3 , F. Niebergall 11 , C. Niebuhr 10 , O. Nix 14 , 

G. Nowak 6 , J.E. Olsson 10 , D. Ozerov 23 , V. Panassik 8 , C. Pascaud 26 , G.D. Patel 18 , M. Peez 22 , 

E. Perez 9 , A. Petrukhin 35 , J.P Phillips 18 , D. Pitzl 10 , R. Poschl 26 , 1. Potachnikova 12 , B. Povh 12 , 
G. Radel 1 , J. Rauschenberger 11 , P. Reimer 29 , B. Reisert 25 , D. Reyna 10 , C. Risler 25 , E. Rizvi 3 , 
P. Robmann 37 , R. Roosen 4 , A. Rostovtsev 23 , S. Rusakov 24 , K. Rybicki 6 , D.P.C. Sankey 5 , 

S. Schatzel 13 , J. Scheins 1 , F.-P. Schilling 10 , P. Schleper 10 , D. Schmidt 33 , D. Schmidt 10 , 
S. Schmidt 25 , S. Schmitt 10 , M. Schneider 22 , L. Schoeffel 9 , A. Schoning 36 , T. Schorner 25 , 
V. Schroder 10 , H.-C. Schultz-Coulon 7 , C. Schwanenberger 10 , K. Sedlak 29 , F. Sefkow 37 , 
V. Shekelyan 25 , 1. Sheviakov 24 , L.N. Shtarkov 24 , Y. Sirois 27 , T. Sloan 17 , P. Smirnov 24 , 
Y. Soloviev 24 , D. South 21 , V. Spaskov 8 , A. Specka 27 , H. Spitzer 11 , R. Stamen 7 , B. Stella 31 , 



1 



J. Stiewe 14 , 1. Strauch 10 , U. Straumann 37 , M. Swart 14 , S. Tchetchelnitski 23 , G. Thompson 19 , 

P.D. Thompson 3 , N. Tobien 10 , F. Tomasz 14 , D. Traynor 19 , P. Truol 37 , G. Tsipolitis 10 ' 38 , 

I. Tsurin 35 , J. Turnau 6 , J.E. Turney 19 , E. Tzamariudaki 25 , S. Udluft 25 , A. Uraev 23 , M. Urban 37 , 

A. Usik 24 , S. Valkar 30 , A. Valkarova 30 , C. Vallee 22 , P. Van Mechelen 4 , S. Vassiliev 8 , 

Y. Vazdik 24 , A. Vest 1 , A. Vichnevski 8 , K. Wacker 7 , J. Wagner 10 , R. Wallny 37 , B. Waugh 21 , 

G. Weber 11 , D. Wegener 7 , C. Werner 13 , N. Werner 37 , M. Wessels 1 , G. White 17 , S. Wiesand 33 , 

T Wilksen 10 , M. Winde 35 , G.-G. Winter 10 , Ch. Wissing 7 , M. Wobisch 10 , E.-E. Woehrling 3 , 

E. Wunsch 10 , A.C. Wyatt 21 , J. Zacek 30 , J. Zalesak 30 , Z. Zhang 26 , A. Zhokin 23 , F. Zomer 26 , and 

M. zur Nedden 10 



1 1. Physikalisches Institut der RWTH, Aachen, Germany 11 

2 Physikalisches Institut der RWTH, Aachen, Germany 11 

3 School of Physics and Space Research, University of Birmingham, Birmingham, UK b 

4 Inter- University Institute for High Energies ULB-VUB, Brussels; Universiteit Antwerpen 
(UIA), Antwerpen; Belgium 

5 Rutherford Appleton Laboratory, Chilton, Didcot, UK b 

6 Institute for Nuclear Physics, Cracow, Poland 8 

7 Institut fur Physik, Universitdt Dortmund, Dortmund, Germany 11 

8 Joint Institute for Nuclear Research, Dubna, Russia 

9 CEA, DSM/DAPNIA, CE-Saclay, Gif-sur-Yvette, France 

10 DESY, Hamburg, Germany 

11 Institut fur Exp erimentalphysik, Universitdt Hamburg, Hamburg, Germany 

12 Max-Planck-Institut fiir Kernphysik, Heidelberg, Germany 

13 Physikalisches Institut, Universiteit Heidelberg, Heidelberg, Germany a 

14 Kirchhoff-Institut fur Physik, Universitdt Heidelberg, Heidelberg, Germany a 

15 Institut fur experimentelle und Angewandte Physik, Universitdt Kiel, Kiel, Germany 

16 Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak Republic 6 ^ 

17 School of Physics and Chemistry, University of Lancaster, Lancaster, UK b 

18 Department of Physics, University of Liverpool, Liverpool, UK b 

19 Queen Mary and Westfield College, London, UK b 

20 Physics Department, University of Lund, Lund, Sweden 9 

21 Physics Department, University of Manchester, Manchester, UK 11 

22 CPPM, CNRS/IN2P3 - Univ Mediterranee, Marseille - France 

23 Institute for Theoretical and Experimental Physics, Moscow, Russia 1 

24 Lebedev Physical Institute, Moscow, Russia 6 

25 Max-Planck-Institut fur Physik, Munchen, Germany 

26 LAL, Universite de Paris-Sud, IN2P3-CNRS, Orsay, France 

27 LPNHE, Ecole Poly technique, IN2P3-CNRS, Palaiseau, France 

28 LPNHE, Universites Paris VI and VII, IN2P3-CNRS, Paris, France 

29 Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic 6 ' 1 

30 Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic 6 ' 1 

31 Dipartimento di Fisica Universita di Roma Tre and IN FN Roma 3, Roma, Italy 

32 Paul Scherrer Institut, Villigen, Switzerland 

33 Fachbereich Physik, Bergische Universitdt Gesamthochschule Wuppertal, Wuppertal, 
Germany 

34 Yerevan Physics Institute, Yerevan, Armenia 



2 



35 DESY, Zeuthen, Germany 

36 InstitutfUr Teilchenphysik, ETH, Zurich, Switzerland^ 

37 Physik-Institut der Universitdt Zurich, Zurich, Switzerland^ 

38 Also at Physics Department, National Technical University, Zografou Campus, GR-15773 
Athens, Greece 

39 Also at Rechenzentrum, Bergische Universitdt Gesamthochschule Wuppertal, Germany 

40 Also at InstitutfUr Experimentelle Kernphysik, Universitdt Karlsruhe, Karlsruhe, Germany 

41 Also at Dept. Fis. Ap. CINVESTAV, Merida, Yucatan, Mexico k 

42 Also at University of P.J. Safdrik, Kosice, Slovak Republic 

43 Also at CERN, Geneva, Switzerland 

44 Also at Dept. Fis. CINVESTAV, Mexico City, Mexico k 

a Supported by the Bundesministerium fiir Bildung und Forschung, FRG, under contract 
numbers 05 HI 1GUA /I, 05 HI 1PAA /I, 05 HI 1PAB /9, 05 HI 1PEA /6, 05 HI 1VHA /7 and 
05 HI 1VHB/5 

b Supported by the UK Particle Physics and Astronomy Research Council, and formerly by the 

UK Science and Engineering Research Council 

c Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT 

d Partially Supported by the Polish State Committee for Scientific Research, grant no. 
2P0310318 and SPUB/DESY/P03/DZ-1/99 and by the German Bundesministerium fUr Bildung 
und Forschung 

e Supported by the Deutsche Forschungsgemeinschaft 

f Supported by VEGA SR grant no. 2/1169/2001 

9 Supported by the Swedish Natural Science Research Council 

1 Supported by the Ministry of Education of the Czech Republic under the projects 

ING0-LA1 16/2000 and LN00A006, by GAUK grant no 173/2000 

i Supported by the Swiss National Science Foundation 

k Supported by CONACyT 

1 Partially Supported by Russian Foundation for Basic Research, grant no. 00-15-96584 



3 



1 Introduction 



The Standard Model of particle physics is known to contain certain processes which violate the 
conservation of baryon and lepton number (B + L) in the case of electroweak interactions and 
chirality in the case of strong interactions [O. Such anomalous processes cannot be described 
by standard perturbation theory. They are induced by instantons [Jj], 0] . In the strong sector, 
described by quantum chromodynamics (QCD), instantons are non-perturbative fluctuations 
of the gluon field. They represent tunnelling transitions between topologically non-equivalent 
vacua. Deep-inelastic scattering (DIS) offers a unique opportunity [^]] to discover a class of 
hard processes induced by small QCD instantons. The rate is calculable]] within "instanton- 
perturbation theory" and is found to be sizeable [|5|-f7|]. Moreover, the instanton-induced final 
state exhibits a characteristic signature [f^|8|-|IT|]. Detailed reviews are given in Refs. [|T2|,|13|]. 



An experimental observation of instanton-induced processes would constitute a discov- 
ery of a basic and novel non-perturbative QCD effect at high energies. The theory and phe- 
nomenology for the production of instanton-induced processes at HERA in positron proton 
collisions at a centre of mass energy of 300 GeV has recently been worked out by Ringwald 
and Schrempp J|,|5|-§]. The size of the predicted cross-section is large enough to make an ex- 
perimental observation possible. The expected signal rate is, however, still small compared to 
that from the standard DIS process. The suppression of the standard DIS background is there- 
fore the key issue in this analysis. QCD instanton-induced processes can be discriminated from 
standard DIS by their characteristic hadronic final state signature, consisting of a large number 
of hadrons at high transverse energy emerging from a "fire-ball"-like topology in the instanton 
rest system [J3],[8|,|9p. Derived from simulations studies characteristic observables are exploited to 
identify a phase space region where a difference between data and the standard DIS simulations 
would indicate a contribution from instanton-induced processes. 

Upper cross-section limits on instanton-induced processes based on standard hadronic final 
state observables measured with the HI detector have already been derived in Refs. [inHT^]. 



Here, we present for the first time a dedicated experimental search for instanton-induced pro- 
cesses in high energy particle collisions. 



2 Phenomenology of QCD Instanton-induced Processes in 
DIS 

Instanton (7) processes in DIS at HERA are discussed within the framework of the work of 
Ringwald and Schrempp [T3],|5|-|9p. These processes dominantly occur in a photon gluon (7*7) 
fusion process as sketched in Fig. [I]. The characteristic /-event signatures result from the fol- 
lowing basic reaction: 

7* + <7^$>fc + g*)+ n g g, (1) 
n f 

'For an exploratory calculation of the instanton-induced contribution to the gluon-induced part of the total DIS 
cross-section at large values of the Bjorken scaling variable x > 0.3, see Ref. [Eh. 
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where g, qp (qp) denotes gluons, right-handed quarks (anti-quarks), rif is the number of quark 
flavours and n g is the number of gluons produced. Right-handed quarks are produced in I- 
induced processes, left-handed quarks are produced in anti-instanton (J) processes. The final 
state induced by instantons or anti-instantons can only be distinguished by the chirality of the 
quarks. Experimental signatures sensitive to instanton-induced chirality violation are not ex- 
ploited in this analysis. Both /-processes and /-processes enter in the calculation of the total 
cross-section. 




Figure 1 : Kinematic variables of the dominant /-induced process in DIS. The virtual photon (4- 
momentum 7 = e — e'), emitted by the incoming positron e, fuses with a gluon (4-momentum 
g) radiated from the proton (4-momentum P). The gluon carries a fraction £ of the longitudinal 
proton momentum. The virtual quark entering the instanton subprocess has 4-momentum q', 
while the outgoing quark (= current quark) from the photon splitting process has q" . Wj is the 
invariant mass of the quark gluon (q'g) system and W is the invariant mass of the total hadronic 
system (the 7P system), s is the invariant mass squared of the 75 system. 

As shown in Fig. [I], a photon splits into a quark anti-quark pair in the background of an 

instanton or an anti-instanton field. The so-called /-subprocess q' + g ^4 X is produced by 
the quark or the anti-quark fusing with a gluon g from the proton. The respective partonic 
final state includes 2 n/ — 1 light quarks and anti-quarks. Therefore, together with the current 
quark (q"), in every /-event, quark anti-quark pairs of each of the nf(— 3) (light) flavours are 
simultaneously produced^. In addition, a mean number of (n g ) ~ 0(l/a s ) ~ 3 gluons is 
expected to be emitted in the /-subprocess. 

The quarks and gluons emerging from the /-subprocess are isotropically distributed in the 
/-rest system defined by q' + g = 0. One expects therefore a pseudo-rapidity^] (77) region with a 

2 In principle, also heavy flavours contribute whenever very small instantons are probed. In general, however, 
the quarks must appear approximately massless on the scale of the dominant effective /-size p e g (Q' 2 , x'), i.e. 
p e s m q <§C 1, where m q is the quark mass. In the HERA kinematic region, the rate is dominated by p e g « 0.35 fm 
such that only up, down and strange quarks appear massless (rif = 3). The contribution of charm and bottom 
quarks to the cross-section is likely to be small. It was checked that the predicted final state signature does not 
change significantly if heavy quarks are included in the simulation. 

3 The pseudo-rapidity of a particle is defined as r\ = — lntan(0/2), where 8 is the polar angle with respect to 
the proton direction defining the +z-axis. 
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width of typically 2 units in r\. This region is densely populated with particles of relatively high 
transverse momentum which are homogeneously distributed in azimuth in the J-rest frame. 
Apart from this pseudo-rapidity band, the hadronic final state exhibits a current jet emerging 
from the outgoing current quark q" . The large number of partons emitted in the /-process leads 
to a high multiplicity of charged and neutral particles in every event. 

The actual number of produced hadrons and their energies crucially depends on the centre 
of mass energy Wj available in the /-system, which in turn can be written (see Fig. [T]) in terms 
of the variables Q' 2 and x' describing the kinematics of the /-subprocess. These variables are 
defined in analogy to the Bjorken scaling variables x and Q 2 . A knowledge of the distributions 
of these variables is indispensable for the correct prediction of the hadronic final state. These 
distributions can be calculated within /-perturbation theory []5],p for large Q' 2 and x'. 

The total /-production cross-section at HERA, <J^ RA , is essentially determined by the 
cross-section of the /-subprocess q' + g ^> X denoted by a^P. The latter can be calculated 
by integrating^ over the / (/)-size p (p) and the // distance 4- vector 



Q' 2 ) 



/POO poo 
d 4 R e i( g+q ')-R dp dp e~ {p+ ~ p)Q ' D(p)D(p) ... e 
Jo Jo 



4tt 



ats(fj. r ) 



n(R 2 /pp,p/ P ) 



(2) 



where several parts of the integrand have been omitted. D(p) (D(p)) is the /-size (/-size) 
distribution that is calculable within /-perturbation theory [HJ] for a s (p r ) In (pp r ) <C 1 with 
a s (p r ) being the strong coupling taken at the renormalisation scale p r and N c = 3 for QCD [[]], 

urn- 



Dip) = d 



2tt 



a s (p r 



2N C 



2tt 
Qs(Mr) 



(p r p) 



ln f +0(a a ) 



(3) 



where / is a known scheme dependent constant. 

The function {l(R 2 / pp, pj p), with —1 < £l(R 2 / pp, pj p) < 0, describes the //-interaction 
associated with a resummation of final state gauge bosons. It is calculable in /-perturbation 
theory, formally for R 2 / pp ^> 1, and may attenuate to some extent the exponent —2n/a s of 
the exponential in equation (Q) that is typical for tunnelling transitions. For a general SU(Nc) 
gauge theory with coupling a, equations (Q) and (Q) give the qualitative behaviour for the /- 
cross-section: 



a 



(I) 
q'g 



2tt 



a 



-£(i+n) 



(4) 



Thus, in the absence of final state gauge boson effects (i.e. Q = 0), with typical values of 
a s ~ 0.4 at HERA and the weak gauge coupling a w 0.033, equation (Q) illustrates the strong 
suppression of electroweak instanton effects which is absent in QCD: 



2tt 



12 



5 > 



2tt 



Qf,, 



7-10 



-148 



(5) 



4 For simplicity, the additional integration over the relative // colour orientation has already been performed in 
equation Both instanton and anti-instanton degrees of freedom enter in the cross-section formula, since it is 
obtained from taking the modulus squared of the amplitude depicted in Fig. ||. The complete formula and more 
details can, for instance, be found in Ref. ^ which contains as well the explicit physical interpretation of the 
variables p and R. 
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In this picture the tiny instanton-induced electroweak B + L violation will only be observable, if 
the final state emission of a huge number of VT-bosons counteracts the exponential suppression 
[p"9|], i.e. (1 + 0) ~ 0. In QCD, however, final state gluons are expected to only provide a 



moderate numerical correction of the rate. Correspondingly, the predictions of the /-induced 
rate in QCD depends much less on the resummation of final state gauge bosons. 

According to equation (Q) the /-size distribution follows a power law: 

D(p) ~ p6-2/3n /+ OK) (6) 

and the integral over p (p) generally diverges for large p (p). However, in the DIS regime 
the exponential factor e~( p+ ^ Q appearing in equation (0) ensures the convergence of the inte- 
gral. For large enough Q' effectively only small size instantons contribute to the cross-section. 
Therefore the /-cross-section is calculable in DIS [§]. 

The p and the Rj p distributions can be calculated using quenched (for nj = 0), non- 
perturbative lattice simulations of the QCD vacuum. They will be discussed and compared 
with perturbative predictions in section [7[ By confronting /-perturbation theory with these lat- 
tice simulations, limits on the validity of /-perturbation theory have been derived [||[7],[Kj]. The 
calculations agree for p < 0.35 fm and Rj p > 1.05, which can be translated into limits on x' and 
Q' 2 , i.e. Q'/Aj^ £ 30.8 and x' > 0.35 [§,[Kjp where is the QCD scale in the ~MS scheme 
for rif flavours. 

For the region, Q' 2 > (30.8 A§L) 2 = 64 GeV 2 , x' > 0.35, x > 10~ 3 and 0.1 < y < 0.9, 
the /-cross-section at HERA has been estimated to be cthera ~ 126 pb [Q. This result has 



recently been updated Ql 1|,|12|,[20|] by using the 1998 world average of the strong coupling [ |21| ] 
to cr^, RA = 89+]J pb. The quoted errors for the /-induced cross-section <J^ RA only contain 
the uncertainty obtained from varying the strong coupling. The change in A^L leads to a 
change of the minimal required Q' 2 to Q'^ n = 113 GeV 2 . This cross-section has been derived 
for three-flavours, corresponding to A^L = 346^29 MeV. An additional cut, Q 2 > Q'^, is 
advocated [§, [10|, [TTJ] to reduce remaining theoretical uncertainties connected with non-planar 
diagrams. 

In this domain the cross-section is o"hera = ^-7°5 P D - The calculation is based on a two 
loop renormalisation group invariant expression of the /-density -D(p) and thus does not depend 
much on the chosen renormalisation scale. In the kinematic domain in which this pioneering 
analysis is performed, i.e. the polar angle of the scattered positron 9 e > 156°, 0.1 < y < 0.6, 
x > 10~ 3 and 10 < Q 2 < 100 GeV 2 , the cross-section calculated with QCDINSQ is cr{^ RA = 
43 pb. 

Even though these predictions have not yet reached the same quantitative level of precision 
as current standard perturbative QCD calculations, the cross-section is large enough to motivate 
dedicated searches for /-processes at HERA. 

5 In this result theoretical uncertainties connected with non-planar diagrams are not taken into account. How- 
ever, the observables used in this analysis to calculate cross-section limits seem to be rather insensitive to these 
uncertainties 
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3 The HI Detector at HERA 



A detailed description of the HI detector can be found elsewhere [|22|]. Here we briefly introduce 
the detector components most relevant for this analysis: the liquid argon (LAr) calorimeter, the 
backward lead-fibre calorimeter (SpaCal) and the tracking chamber system. 



The hadronic energy flow is mainly measured by the LAr calorimeter [|23|] extending over 
the polar angle range 4° < 9 < 154° with full azimuthal coverage. It consists of an elec- 
tromagnetic section (20 — 30 radiation lengths) with lead absorber and a hadronic section 
with steel absorber. The total depth of both calorimeters varies between 4.5 and 8 interac- 
tion lengths. Test beam measurements of the LAr calorimeter modules show an energy reso- 



lution of o E /E « 0.12/ y/E [ GeV ] © 1% for electromagnetic showers [E^] and of a E /E 



0.50/ y/E [ GeV ] © 2% for charged pions 

The backward lead-fibre calorimeter SpaCal [ |2~6l ] covers the polar angle range 153° < 
9 < 177°. In the electromagnetic section, with a depth of 28 radiation lengths, the posi- 
tion and the energy of electrons are measured. The electron energy resolution is o E /E 



0.075/ y/E [GeV] © 2.5%. In total, the SpaCal has two interaction lengths which provide 
additional measurements for hadrons. 

The calorimeters are surrounded by a superconducting solenoid providing a uniform mag- 
netic field of 1.15 T parallel to the beam axis for momentum measurement of charged particles. 
These are measured in two concentric jet drift chamber modules (CJC), covering the polar angle 
range 15° < 9 < 165° [|2"7|]. A backward drift chamber (BDC) aids identification of positrons 
scattered into the SpaCal calorimeter. 

The luminosity is measured using the elastic Bethe-Heitler process ep — > ep'j. The final 
state positron and photon are detected in calorimeters situated close to the beam pipe at distances 
of 33 m and 103 m from the interaction point in the positron beam direction. 



4 Simulation of Standard DIS and /-Processes 



Detailed simulation of the HI detector response to hadronic final states have been performed 
for QCD models of the standard DIS processes and for QCD /-induced scattering processes. 



4.1 Simulation of Standard DIS 



The RAPGAP Monte Carlo [|2|] incorporates the 0(a s ) QCD matrix element and models 
higher order parton emissions to all orders in a s using the concept of parton showers p3|] 



based on the leading logarithm DGLAP equations Q30Q, where QCD radiation can occur before 
and after the hard subprocess. The formation of hadrons is performed using the LUND string 
model [jnjl implemented in JETSET [|32p. This QCD Monte Carlo is called "MEPS". 



An alternative treatment of the perturbative phase is implemented in ARIADNE [B3J], where 



gluon emissions are simulated using the colour dipole model (CDM) [34] by assuming a chain 
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of independently radiating dipoles spanned by colour connected partons. The first emission in 



the cascade is modified to reproduce the matrix element to first order in a s Q35|]. The hadroni 



sation is performed using JETSET This QCD Monte Carlo is called "CDM" in the following. 
The CDM and the MEPS Monte Carlo simulations are both interfaced to the program HER- 



ACLES pq ] to include 0(a) electroweak corrections to the lepton vertex, where a is the elec- 



tromagnetic coupling. 

The Monte Carlo DIS samples have been generated using the CTEQ4 [TJBp parton density 



functions and have been reweighted using a parametrisation extracted from the recent HI mea- 



surement of the proton structure function [ ]39(1 . 



The MEPS and the CDM models have been comprehensively compared to a variety of 
hadronic final state data, and an attempt has been made to optimise the free model parameters 
pO|,|4T|]. No single parameter set was found which describes all studied distributions well. 
Moreover, the Monte Carlo models make different predictions. Therefore and in view of the 
involved approximations it is questionable to what extent the currently available QCD models 
can describe the standard DIS hadronic final state in particular in the tails of distributions. In 
this analysis the Monte Carlo models have been used with their default parameter values. 



4.2 Simulation of QCD Instanton-induced Processes 

QCDINS [|TT|,[42|] is a Monte Carlo package to simulate QCD /-induced scattering processes in 
DIS. It acts as a hard process generator embedded in the HERWIG [ ft3| ] program. The hard pro- 
cess is treated according to the physics assumptions explained in section 0. Apart from the Q 2 
cut, the default parameters of the QCDINS 2.0 version were used, i.e. x' > 0.35, Q' 2 > 113 
GeV 2 and the number of flavours is set to nj = 3. The CTEQ4 parton density functions 
have been employed. After assembling the hard J-subprocess, further QCD emissions are sim- 
ulated in the leading-logarithm approximation. The coherent branching algorithm implemented 
in HERWIG is used. The transition from partons to the observable hadrons is performed with 
the cluster fragmentation model (fflp. 

The hadronic final state topology is mainly influenced by the energy available for the hard 
J-subprocess. It has been explicitly checked that the conclusions drawn from this analysis are 
unchanged when the LUND string model is used instead of the cluster fragmentation model. 
This has also been observed in Ref. [Q], where in addition the effect of changing free model 
parameters in the hadronisation models has been studied. For what follows it is assumed that 
the commonly used hadronisation models are also applicable to describe the fragmentation of 
a large number of (9(10) partons produced by the I-process in a narrow pseudo-rapidity region 
with high transverse energy. 



5 Event Selection and Search Strategy 
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5.1 Inclusive DIS Event Selection 



The data used in this analysis were collected in the years 1996 and 1997 with the HI detector 
at the electron proton collider HERA. During this time HERA collided positrons at an energy 
of E e = 27.5 GeV with protons at an energy of E p = 820 GeV. The accumulated data sample 
corresponds to an integrated luminosity of 21.1 pb -1 . 

The scattered positron is identified as the electromagnetic energy deposition with the highest 
energy. For this pioneering analysis we restrict the electron energy and angle measurement to 
the SpaCal calorimeter. The electron is required to lie well within the calorimeter and trigger ac- 
ceptance of polar angles between 156° and 176°. A minimal positron energy of E' e > 10 GeV is 
required. The events are triggered by demanding a localised energy deposition in the SpaCal 
together with loose track requirements in the multi-wire proportional and the drift chambers. 

Furthermore the longitudinal momentum balance is required to lie within 35 GeV < 
p z ) < 70 GeV, where the sum runs over the scattered electron and all objects belonging to the 
hadronic final state. The hadronic final state objects are reconstructed from the calorimetric 
energy depositions in the LAr and the SpaCal calorimeters and from low momentum tracks 



(0.15 < p t < 2 GeV) in the central jet chamber according to the procedure described in [|4q]. 
The position of the z coordinate of the reconstructed event vertex must be within ±30 cm of the 
nominal interaction point. 

The photon virtuality Q 2 and the Bjorken scaling variable x are reconstructed from the 
scattered positron. The events are selected to cover the phase space region defined by 9 e > 156°, 
0.1 < y < 0.6, x > 10~ 3 and 10 < Q 2 < 100 GeV 2 , where 6 e is the polar angle of the scattered 
positron. 

The selected DIS data sample consists of about 375000 events. The simulated events repro- 
duce well the shape and the absolute normalisation of the distributions of the energy and angle 
of the scattered positron as well as the kinematic variables x, Q 2 and y. The contamination with 
events due to hadrons misidentified as positrons produced in collisions of high energetic protons 
with quasi-real photons is below 2%. This was estimated using the Monte Carlo simulation pro- 
gram PHOJET p7|], which contains the 0(a 3 ) matrix elements for direct and resolved photon 
processes, parton showers and a phenomenological description of soft interactions. 



5.2 Definition of the Discriminating Observables 

The observables used to discriminate the /-induced contribution from the standard DIS process 
are based on the hadronic final state. Only hadronic final state objects as defined in section 



5TT| within —1.4 < < 2 are considered. Charged particles with transverse momenta of 
p t > 0.15 GeV are selected within 20° < 9 < 155°. Here, both r]^ and p t are measured in the 
laboratory frame. 

All hadronic final state objects are boosted to the hadronic centre-of-mass frameQ. Jets 
are defined by the cone algorithm [|46|, |47|] with a cone radius of R = a/ At] 2 + A0 2 = 0.5. 



6 The hadronic centre-of-mass frame is defined by 7 + P = 0, where 7 (P) is the 3-momentum of the exchanged 
photon (proton). 
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The jet with the highest transverse energy (E t j et ) is used to estimate the 4-momentum q of 
the current quark (see Fig. [[]). Q' 2 can be reconstructed from the particles associated with the 
current jet and the photon 4-momentum, which is obtained using the measured momentum of 
the scattered positron. The cone size was chosen to optimise the resolution of the reconstruction 



of Q' (see also Ref. [fTS]). The Q' resolution is about 20 — 30%. However, the distribution 



of the true over the reconstructed value exhibits large tails, since in about 30% of the cases the 
wrong jet is identified as the current jet. Due to the limited accuracy of the Q' 2 reconstruction, 
the reconstructed Q' 2 cannot be used to experimentally control the "true" Q' 2 region of the I- 
processes, but can nevertheless be exploited to discriminate /-processes from the standard DIS 
background. The reconstructed Q' 2 is called Q' V 2 C in what follows. More information on the 
Q' 2 reconstruction can be found in [§, ^8[^]. The reconstruction of the variable x' is more 
difficult[| an d on ly possible for /-events exhibiting the characteristic /-topology [|49|]. 

The hadronic final state objects belonging to the current jet are not used in the definition of 
the following observables. A band in pseudo-rapidity with a width of 1.1 units in r] is defined 
around the centre of gravity fj = ^ E T i]/ E T ) of the transverse energy (E T ) distribution of 
the hadronic final state objects (see Ref. p9[ ] for details). This pseudo-rapidity band is called 
the /-band in the following. The number of charged particles in the /-band measured as tracks 
in the detector is counted (n B ) and the total scalar transverse energy of all hadronic final state 
objects in the /-band is measured (E^b)- 

All hadronic final state objects in the /-band are boosted to an approximate /-rest frame 
defined by q' + (£)P = 0, where (£) = 0.076 is the average value expected by the QCDINS 
Monte Carlo simulation (see Fig. [I] for definition). In this system the sphericity (Splie) is 
calculated^. For spherical events Sph B is close to 1, while for pencil-like events Sph B is 0. 
Furthermore, the axes i min and i max are found for which in the /-rest system the summed pro- 
jections of the 3-momenta of all hadronic final state objects in the /-band are minimal or max- 
imal [§]. The relative difference between E in = Y,h Wh ' «max| and E out = J2h \Ph ' w| is 
called Ab = (E in — E out )/E in . This quantity is a measure of the transverse energy weighted 
azimuthal isotropy of an event. For isotropic events A B is small while for pencil-like events 
A B is large. 

Three observables are used to enhance the fraction of /-events in the inclusive data sample: 
the charged particle multiplicity in the /-band (n B ), the sphericity of all final state objects in 
the /-band calculated in the approximate /-rest frame (Sphs) and the reconstructed Q{ 2 C . Three 
other observables, i.e. E tyJet , E tjB and A B , contain additional information and will be used for 
further checks. 



7 The data used for the present analysis do not incorporate a cut on x' . As noted in Refs. [0,|l^], this pre- 
sumably does not prevent a qualitative comparison with the Ringwald-Schrempp predictions in the fiducial x' and 
Q' 2 region for the following reason. The lattice data for the //-distance R distribution (from which the min- 
imal theoretical x'-cvX was deduced [7],[k]]) exhibit a rapid suppression of /-effects for small //-separation, 
corresponding to x' < 0.35. Therefore, /-contributions to the data from this x' region outside the validity of 
/-perturbation theory can probably be neglected. 

8 The sphericity is defined as SPH = (3/2)(A2 + A3) where A2 and A3 are the smallest of the three eigenvalues 
of the diagonalised sphericity tensor defined by S a @ = (X^Pfff )/ Si b«| 2 ' where a and /3 corresponds to the 
x, y and z components of the considered particle momenta pi. 
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5.3 Comparison of Data to Standard QCD Predictions in Inclusive DIS 



The distributions of the observables Splie and Q.' ec for data, for two standard DIS QCD 
models and for the /-process are shown in Figs. ^ and [3|. The gross features of the data are 
reasonably well described by both Monte Carlo simulations. The CDM model is able to describe 
the data within 10% except at very low and very large sphericity values where a difference of 
20% is observed. The MEPS Monte Carlo reproduces the data within 10 — 15%. However, at 
large n B deviations up to 30% are found. 

The three other observables, i.e. E tJet , E tB and A B , used as control distributions are shown 
in Figs. and J3L The A B distribution is fairly well described (within 10 — 20%) by both standard 
DIS simulations. For E t; j et the two standard DIS Monte Carlos simulations behave differently. 
MEPS describes the data within 10% for E t:Jet < 2.5 GeV, but is lower by about 20% for 
E t j et > 2.5 GeV. In the tail of this distribution (for E t „ Jet > 10 GeV) the data are well 
reproduced up to the largest accessible values. The CDM model describes the data within 
5 — 10% for E tJet < 5 — 10 GeV, but the predictions progressively grow to be above the data 
at higher values. A harder tail than found in the data is also seen in the E tjB distribution. CDM 
overshoots the data by 50% at large E tjB values. The MEPS simulation gives a much better 
description of this observable (within 20%). The hard transverse energy tail produced by CDM 
has also been observed in two-jet and three-jet production in DIS at HERA pT], pT||52[ ]. 

The instanton prediction is shown as a dotted line in Fig. ^ and Fig. |3|. For visibility in Fig. |2| 
it is scaled up by a factor 500. The expected /-signal is about two to three orders of magnitude 
smaller than the standard DIS background. Therefore cuts are needed to enhance the signal to 
background ratio. 



6 Search for Instanton-Induced Events 

Two methods are employed to increase the sensitivity to /-processes: a combinatorial cut 
method described in section lO and a multivariate discrimination technique based on a range 



search method described in section pJ\ While the main advantage of the combinatorial cut 
method is its simplicity, the multivariate discrimination technique allows the transition from a 
background dominated to an /-enriched region to be investigated with one single observable. 



6.1 Combinatorial Cut Based Method 

The strategy to reduce the standard DIS background is based on the observables n B , Q{. ec and 
Sph B . These observables have been chosen, since they provide the best signal to background 
separation. Moreover, in the case of /-induced processes the shape of their distributions is 
expected not to be much affected by contributions from non-planar diagrams [ ]T0| ] (see also 
section 0). 

To find the optimal combinations of cut values the three observables are investigated using 
simulations of the standard DIS background and of the /-signal. Amongst the studied cut com- 
binations with an instanton efficiency e/>10% the one is chosen with the best separation power, 
defined by S = ei/e s nis, where e S Dis is the fraction of remaining standard DIS background. 
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Combinatorial cut method (section |6.1 ) 




N 


(■sDIS 






Data 
CDM 
MEPS 


484 

443+1 

304j£ 


0.118% 
0.081% 


86 
125 


47 pb 
109 pb 


Multivariate method (section ^ 


2) 




N 


e sDIS 


£// e sD/S 




Data 
CDM 
MEPS 


410 

354+*° 


0.095% 
0.080% 


106 
126 


55 pb 
80 pb 



Table 1: Number of events observed in the data and expected from the CDM and MEPS simu- 
lation after optimising the /-signal to background ratio. The quoted error contains the full sta- 
tistical and systematic uncertainty added in quadrature. The efficiency of the standard DIS sim- 
ulation (e S Dis), the separation power (S = £i/e S Dis) an d the exclusion limit at 95% confidence 
level on the production cross-section of /-induced processes simulated in the fiducial region 
x' > 0.35 and Q' 2 > 113 GeV 2 in the kinematic region defined by x > 10~ 3 , 0.1 < y < 0.6, 
Q 2 < 100 GeV 2 and# e > 156° are also given. The /-cross-section calculated by /-perturbation 
theory is 43 pb. 



The best separation power S is obtained for 95 < Q^ c < 200 GeV 2 , n B > 11 and Sph B > 
0.4. For this cut combination S = 125 (S = 86) is found for the MEPS (CDM) simulation. At 
an /-efficiency of about 10%, the background has been suppressed by about a factor of 1000. 
With these cuts, 484 events are found in the data, while CDM predicts 44311 and MEPS 304+Ji; 
(see also Tab. [I]). The quoted errors on the expected event numbers include the statistical and 
the experimental systematic uncertainties. 

The uncertainties arise from the following sources: a 4% uncertainty in the hadronic en- 
ergy scale of the LAr calorimeter, 1% for the electromagnetic energy scale and 7% for the 
hadronic energy scale measured in the SpaCal calorimeter, 3% for the measurement of the track 
momentum, 2 mrad for the polar and azimuthal angle of the track, 2% (5%) for tracks with a 
momentum above (below) 250 MeV for inefficiencies in the track reconstruction and 2 mrad 
for the polar angle of the scattered positron. These uncertainties have been propagated into 
the overall systematic error. An absolute normalisation error of 1.5% for the accuracy of the 
luminosity determination and a 3% uncertainty for the reweighting of the parton densities used 
in the Monte Carlo simulation have been included. The main contributions to the systematic 
uncertainties on the number of expected events are seen to arise from the track efficiency (3.5%) 
and the momentum measurement (3%). 

More events are found in the data than expected by either one of the background Monte 
Carlo simulations. MEPS suggests a clear excess in the data. In contrast to the findings before 
placing the cuts CDM is compatible with the data within errors. The predictions of the two 
standard DIS Monte Carlo simulations largely disagree with each other. This indicates that the 
background estimation is subject to large uncertainties as it is expected in this extreme region 
of phase space. The distribution of the discriminating observables after the cuts is shown in 
Fig. |] in comparison with the standard DIS background expectation from the MEPS and the 
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CDM Monte Carlo models. Fig. [| also shows the expectation for the /-contribution to the DIS 
process as modelled by QCDINS. In particular when compared to the MEPS model, the shape 
of the observed excess in the data is qualitatively compatible with an /-signal for n B , Q^ c and 
Splie- However, it tends to lie towards low E t) s values in contrast to the /-contribution. It has 
been noted in [ |T0"| ] that the E^b (as the E t j et ) distribution expected for the /-signal are most 
sensitive to theoretical uncertainties. A shift towards lower values is well possible within these 
uncertainties. 



6.2 Search Based on a Multivariate Discrimination Technique 

To make optimal use of the information contained in the observables separating the /-signal and 
the background a multivariate discrimination technique is used. Events are classified as signal 
or background by estimating their probability density p at each point in the phase space of the 
observables. Monte Carlo simulations are employed to sample these densities. The densities at 
each phase space point can be directly estimated by counting the number of expected signal and 
background events in a surrounding box. The likelihood of an event to be due to a signal can be 
defined by: 

D= & 

p(I) + p(sDIS) 

By cutting on this discriminator D, the background contamination can be minimised at the 
expense of signal efficiency. For each value of D the signal to background ratio as well as the 
signal efficiency can be easily calculated. 

To estimate the p(I) and p(sDIS) in the vicinity of a given event in the multi-dimensional 



phase space, a range search algorithm [ p3| ] based on a binary tree has been used. A more 



detailed description and the properties of this method can be found in [ p4[p5| ]. 

For the discriminator the phase space is spanned by the same observables n B , Q.' ec and 
SpliB , which are used for the combinatorial cut method. The discrimination power of the 
method is demonstrated in Fig. ^a where the shape of the discriminator distribution is shown. 
The simulated background events (solid and dashed line) are mainly concentrated at low D val- 
ues while the simulated /-signal events (dotted line) are peaked towards D — 1. For ej = 10%, 
which corresponds to a cut at D > 0.988, a separation power of S = 126 for MEPS and 
S = 106 for CDM is obtained, respectively. The separation power is slightly improved with 
respect to the combinatorial cut method. 410 events are observed in the data, while 354+26 
(299^3g) are expected for CDM (MEPS) (see also Tab. []]). These results are consistent with the 
ones obtained from the cut-based method. The dominant contributions to the systematic un- 
certainties are attributed to the track efficiency (3.5%), the momentum measurement (3 — 6%) 
and the energy scales in the calorimeters (4 — 5%). The discriminating observables after a cut 
D > 0.988 are shown in Fig. |[ Note that the phase space region selected by the discriminator 
roughly coincides with the results of the combinatorial cut method. The observable distributions 
are suppressed at approximately the same points where the cuts have been placed by the combi- 
natorial cut method. The shape of the excess in the data in the n B , Q^ c and Sph B distributions is 
similar to the shape of the expected /-distributions. The excess in the E t B distribution is largest 
at E t ,B ~ 9 GeV while from /-processes a peak at about E ty B ~ 12 GeV is expected. Again, 



14 



it should be noted that these distributions might be shifted towards lower values when contri- 



butions from non-planar diagrams are taken into account [|10|]. The fact that for the E t B and 
the E t Jet distributions the standard DIS Monte Carlo simulation also disagree with each other 
indicates large uncertainties from the different treatment of higher order QCD processes. 

The multivariate discrimination technique offers furthermore the possibility to compare the 
description of the data by the Monte Carlo simulation in the complete phase space, i.e. from the 
region where no /-contribution is expected (D = 0) to the /-enriched region (D = 1). Fig. |^b 
shows the absolutely normalised discriminator distribution on a logarithmic a;-axis scaled as 
— log 10 (1 — D). The majority of the standard DIS background events are concentrated at the 
lowest D values. Towards larger D values the background falls by three orders of magnitudes. 
The data roughly follow this trend. In the last three bins, a slight excess of data over background 
is observed. According to the QCDINS simulation, in this region 10% of all /-events are con- 
tained. As shown in Fig. |6|d, where the expected /-signal with respect to the data is shown, the 
/-contribution in the event sample is about 20%. 

The description of the data by the background Monte Carlo simulations is illustrated in 
more detail in Fig. fy, where the relative difference of background and data is shown. The 
MEPS Monte Carlo gives an excellent description of the data for D < 0.90 which comprises 
the vast majority of events. Towards larger D values an increasingly large excess of events is 
seen in the data. The largest discrepancy of 60% is found at the largest D value. The excess of 
data as a function of D is qualitatively similar to the increasing ratio of QCDINS events to the 
data (see Fig. f 



From a purely statistical point of view the excess in particular in case of the MEPS back- 
ground estimation is significant. However, the uncertainties in the background estimation are 
largely unknown in this extreme phase space region. This is also reflected in the different 
behaviour of the CDM Monte Carlo simulation. While it agrees with the data in the pure back- 
ground region D < 0.2, it is not able to describe the data at larger values (0.25 < D < 0.9), 
where no significant /-contribution is expected. In this region CDM in contrast to MEPS pre- 
dicts more events than found in the data. It is interesting that both CDM and MEPS fall below 
the data when the separation, i.e. D, is largest. In the case of CDM this observation is, however, 
not significant given the experimental uncertainties. Whether the excess can be explained by 
/-processes or whether it is simply due to a deficiency in the description of the relevant standard 
DIS process remains an open question. Altogether, despite some excess of events in the /-signal 
region, the uncertainties in the background estimation are too large to draw firm conclusions. 



7 Exclusion Limits for Instanton-induced Processes 

Since no significant excess can be claimed upper limits on the QCD instanton cross-section are 
derived. 

The hadronic final state of /-induced events is strongly influenced by the centre-of-mass en- 
ergy squared, Wf = Q' 2 (1 —x' )/x', available for the partons emerging from the /-subprocess. 
The distributions of the final state topology therefore crucially depend on both the minimal cut 
values, above which /-perturbation theory is expected to be valid, and on the assumed x' and 
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Q' distributions, which are motivated by the validity of /-perturbation theory and which are 
only under theoretical control for large enough x' and Q' 2 . 

In this section first upper cross-section limits for instantons produced in the fiducial region 
x' > 0.35 and Q' 2 > 113 GeV 2 are derived. It is assumed that the x' and Q' distributions are 
correctly described by /-perturbation theory. This approach is based on the results described 
in section ||. To be less dependent on theoretical assumptions instantons are independently 
simulated in several bins of approximately constant x' and Q' 2 and the analysis is repeated. In 
this way the uncertainty on the assumed shape of the x' and Q' 2 distributions is minimised and 
the analysis is extended towards lower x' and Q' 2 value. Here, only the reasonable assumption 
is made that the topology of an /-event is for a given Wj correctly modelled by the /-Monte 
Carlo simulation. 

7.1 Exclusion Limits in the /-Fiducial Region 

First an upper limit on the /-cross-section is derived at 95% confidence level (CL) for instanton 
produced in the fiducial region, where the x' and Q' 2 distributions are calculated within the 
Ringwald-Schrempp approach. The number of observed events in the data and of expected 
standard DIS background events and the /-signal efficiency obtained with the combinatorial cut 
method are used. To derive the cross-section limits the method described in Ref. fl56| ] was used. 
For the background estimation and the selection efficiency statistical and systematic errors are 
taken into account by folding Gaussian distributions into the integration of the Poisson law used 
to determine the limit. 

An /-cross-section of 109 pb (47 pb) is excluded when MEPS (CDM) is assumed to pro- 
vide the correct background description. These results together with the ones obtained from the 
multivariate discrimination technique are summarised in Tab. [I]. Both methods lead to similar 
results. The limits are not far from the predicted /-cross-section of about 43 pb. 

As said before, it is questionable whether the CDM and MEPS models are able to adequately 
describe the standard DIS background in this extreme corner of phase space, where only ~ 0.1% 
of the events in the total sample of standard DIS events are expected. To be independent of 
the detailed modelling of the hadronic final state of DIS events, an additional upper limit is 
extracted where the expected standard DIS background is assumed to be zero. Whatever the 
"true" number of standard DIS events in the selected corner of phase space is, the /-signal can 
certainly not be bigger than the number of observed events in the data. An upper limit on the 
/-cross-section derived in this way is therefore the most conservative one, since it only uses the 
expected topology of /-induced events and the number of events observed in the data. Here 
the combinatorial cut method was used. At 95% confidence level a cross-section of 221 pb is 
excluded without relying on the correct modelling of the background. The /-cross-section 
predicted within the Ringwald-Schrempp framework is about a factor of 5 lower. 

7.2 /-Model Independent Exclusion Limits 

To minimise the theoretical input in the extraction of an upper limit on the /-cross-section, 
small ranges of approximately constant Q' 2 and x' are analysed. Events in a 5 x 5 grid with 
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0.2 <x'< 0.45 (grid size 0.05) and with 60 < Q' 2 < 160 GeV 2 (grid size 20 GeV 2 ) have 
been simulated using the QCDINS Monte Carlo program. In this way the effect of the assumed 
shape of the x' and Q' 2 distributions is minimised and the analysis is extended into regions, 
where the x' and Q' 2 distributions cannot be calculated. It is only assumed that the hadronic 
final state for instantons produced at fixed Q' 2 and x', i.e. at fixed Wf, is correctly modelled by 
the /-Monte Carlo simulation. 



For each particular x' and Q' bin, the analysis described in section |6J] is repeated. The 
new best cut combinations ensuring an /-efficiency ej > 10% and the maximal separation 
power ei I £ S Dis are chosen. The data selected for these optimised cuts on Ub, Q^. c an d Splie are 
compared to the standard DIS background simulations. For each case an upper limit on the I- 
cross-section at 95% CL is derived. In Fig. [7] the cross-section exclusion limit obtained with the 
CDM and MEPS background simulations are shown for each x' and Q' 2 region. The difference 
of the results for the two standard DIS models reflects the uncertainty in the DIS prediction. The 
conservative limit independent of the background description, i.e. for the background assumed 
to be zero, is also shown. Depending on the x' and Q' 2 intervals /-cross-sections between 
60 pb and 1000 pb are excluded. 

In Fig. [7] also the /-cross-section as evaluated in the Ringwald-Schrempp framework is 
shown in the region where /-perturbation theory is expected to be valid, i.e. for large x' and 
Q /2 corresponding to small instanton sizes. The upper limits are above the predicted /-cross- 
sections in the fiducial region x' > 0.35 and Q' 2 > 113 GeV 2 . The continuation of the 
fast increase of the /-cross-section towards lower values of x' and Q' 2 as expected from an 
extrapolation of /-perturbation theory is also indicated. 

The prediction of /-perturbation theory can be compared with non-perturbative lattice sim- 
ulations of the QCD vacuum for zero flavours [[7J, to obtain independent information on the 
p and the Rj (p) distributions. This is illustrated^] in Fig. ||a and Fig. ||b. For low p and large 
Rj (p), where the /-perturbation theory is expected to be valid, a power- like rise oc p 6 (see equa- 
tion (§)) is found which is (including the absolute normalisation) in excellent agreement with 
the non-perturbative lattice simulations. For p > 0.35 fm and Rj (p) < 1.05 the /-perturbative 
calculation continues to rise while the lattice data flatten and finally drop towards p — > 1 fm 
and small R/ (p), respectively. In this region the /-perturbation theory is not reliable anymore. 
Figure ||c shows an expanded version of the region where /-perturbation theory starts to deviate 
from the lattice data. 

To compare this observation with the HERA cross-section limits we transform the x' and 
Q' 2 bins into the p and Rj (p) bins by using the effective /-size p e s(x', Q' 2 ), which dominates 
the integration in equation (0) and can be obtained from information accessible in QCDINS. 
To present the cross-section limits we choose one Rj (p) bin which is just outside (0.99 < 
Rj (p) < 1.06) and one bin which is just inside (1.06 < Rj (p) < 1.12) the fiducial region. 

0. 99 < R/(p) < 1.06 corresponds to 0.3 < x' < 0.35 and 1.06 < R/(p) < 1.12 corresponds 
to 0.35 < x' < 0.4. The obtained cross-section limits are shown in two bins of Rj (p) in Fig. ||d 
as a function of p e fj. The QCDINS predictions are shown as lines. For 1.06 < R/(p) < 1.12, 

1. e. within the fiducial region of /-perturbation theory, the cross-section limits exclude the 

9 The notation is as follows: d t^' dp corresponds to D(p) and rf 4 "jl^ = J °° dp / °° dpD(p)D(p) e ~ in 
eq. ||. The variable x denotes the Euclidian space-time coordinates. 
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continuation of the power-like rise of the /-cross-section with p for p > 0.35 fm. For 0.99 < 
R/(p) < 106 the rise has to be attenuated even earlier. This disfavours a continuation of a 
steep rise of the /-cross-section towards large p values. The absence of this rise is in accord 
with the expectation of lattice simulation of the QCD vacuum. 

8 Summary 

QCD instanton-induced processes as modelled in a Monte Carlo simulation implementing the 
prediction of /-perturbation theory have been searched for in deep-inelastic scattering at HERA 
in the kinematic range x > 10~ 3 , 0.1 < y < 0.6, 8 e > 156° and Q 2 < 100 GeV 2 . Three 
observables were used to enhance the sensitivity to instanton events with respect to the standard 
DIS QCD background: the charged particle multiplicity in the instanton rapidity band, the 
reconstructed quark virtuality Q^ c and the sphericity of the hadronic final state objects in the 
instanton rapidity band. Applying either cuts or a multivariate discrimination technique based 
on these three observables, the standard DIS background is suppressed by typically a factor of 
1000, while 10% of the /-events are expected to be kept. In this region 484 events were observed 
and 443l|| and 304± 25 were predicted by two standard DIS QCD background models. 

Using a multivariate discrimination technique the transition from the background dominated 
region to the /-enriched region is investigated. The matrix element plus parton shower model 
describing the data in the background region clearly falls below the data in the signal region. 
With increasing sensitivity to /-processes an increasingly large excess is seen in the data. The 
shape of the excess is qualitatively compatible with the expected /-signal. The colour dipole 
model does not describe the data well in the background dominated region, but is in better agree- 
ment in the signal region. Although the data exceed the expectations where the sensitivity to the 
instanton process is expected to be largest, this effect is not significant given the uncertainties 
in the background estimation. 

The standard DIS Monte Carlo models have known deficiencies and fail to describe vari- 
ous aspects of DIS data in the HERA regime. A better understanding of the formation of the 
hadronic final state in DIS in the bulk and in the extreme end of the phase space relevant for 
instanton searches will be needed to make further progress. 

Based on solely the number of observed events in the data to remain independent of the 
modelling of the standard DIS background, a most conservative upper limit on the instanton 
cross-section of 221 pb is excluded at 95% confidence level. This limit is valid in the fiducial 
region of instanton perturbation theory, i.e. for x' > 0.35 and for Q' 2 > 113 GeV 2 implying 
small /-sizes p and large // distances R, where the /-calculations are expected to be valid. The 
limit is about a factor of five above the cross-section predicted by instanton perturbation theory 
and relies on the /-topology based on the calculated x' and Q' 2 distributions. To be independent 
of this assumption additional upper exclusion limits on the /-cross-section have been derived 
in fixed x' and Q' 2 intervals corresponding to small regions of p and Rj (p). These limits are 
only based on the hadronic final state /-topology for a given x' and Q' 2 and on the number of 
observed events in the data. Depending on the considered kinematic region /-cross-sections 
between 60 and 1000 pb are excluded at 95% confidence limit. The limits cannot exclude the 
predicted /-cross-section in the fiducial region, but exclude a steep /-cross-section rise with 
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decreasing x', i.e. towards large /-sizes, as would be obtained from a naive application of I- 
perturbation theory in this region. The absence of such a steep rise is in accord with lattice 
simulations for zero flavours of the QCD vacuum. 

In summary, this initial experimental study has shown that very large instanton contributions 
at small momentum transfers can be excluded at HERA. To reach sensitivity at the level of the 
predicted instanton-induced cross-section requires improved sophistication in the experimental 
methods and a thorough understanding of the hadronisation process. 
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Figure 2: Distributions of (a) the sphericity in the /-band, Sph B , (b) the reconstructed virtu- 
ality, Q' r l c , (c) the charged particle multiplicity in the /-band, (d) the total transverse energy in 
the /-band, E t>B , (e) the isotropy variable A B and (f) the transverse current jet energy, E t>Jet , 
in the inclusive DIS sample. Data (filled circles), the QCD model background Monte Carlo 
simulations (solid and dashed line) and the QCDINS prediction scaled up by a factor of 500 
(dotted) are shown. 
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Figure 3: Same distributions as in Fig. ^ except for the absence of rescaling the QCDINS 
predictions, but on a logarithmic scale. 
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Figure 4: Distributions of observables after the combinatorial cuts: (a) the sphericity in the 
/-band, Sph B , (b) the reconstructed virtuality of the quark, Q^ c , (c) the charged particle mul- 
tiplicity in the /-band, n B , (d) the total transverse energy in the /-band, E tjB , (e) the isotropy 
variable A B and (f) the transverse current jet energy, E t j et . Data (filled circles), the QCD 
model background Monte Carlo simulations (solid and dashed line) and the QCDINS predic- 
tion (dotted) are shown. 
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Figure 5: Distributions of observables after a cut on the discriminator D > 0.988 to enrich 
/-events: (a) the sphericity in the /-band, Sph B , (b) the reconstructed virtuality of the quark, 
Qrec ( c ) me charged particle multiplicity in the /-band, n B , (d) the total transverse energy in the 
/-band, E t:B , (e) the isotropy variable A B and (f) the transverse current jet energy, E tyJet . Data 
(filled circles), the QCD model background Monte Carlo simulations (solid and dashed line) and 
the QCDINS prediction (dotted) are shown. Data (filled circles), the QCD model background 
Monte Carlo simulation (solid and dashed line) and the QCDINS prediction (dotted) are shown. 
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Figure 6: (a) normalised and (b) event number distribution of the discriminator D and (c) the 
relative difference of data and standard DIS background Monte Carlo simulation. In (d) the ratio 
of the QCDINS prediction to the data is shown. The bands indicate the experimental uncertainty 
of the standard DIS background Monte Carlo. The vertical lines in (c) illustrate the statistical 
error of the data. In (b), (c) and (d) the a;— axis is scaled as — log 10 (1 — D). 
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Figure 7: Upper limit on the cross-section for instanton-induced events as modelled by QCDINS 
for /-events produced in bins of x' and Q' 2 . Regions above the curves are excluded at 95% 
confidence level. The instanton cross-section calculated in the fiducial region x' > 0.35 and 
Q' 2 > 113 GeV 2 (solid line) is also shown and continued towards lower x' (dashed line). The 
limit is valid in the kinematical region of x > 1CT 3 , 0.1 < y < 0.6 and 9 e > 156°. 
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Figure 8: Instanton density distribution on the QCD vacuum as a function of the instanton 
size p (a) and (c) (zoomed) and of the instanton-anti-instanton distance R distribution as a 
function of i?/(p)normalised to the value at R = 1.5 fm (b). Lattice data by the UKQCD 
collaboration [0,^7]] are shown as closed symbols, the prediction of perturbative instanton theory 
as lines [Q]. The dashed vertical lines shows the lower edge of the region where the deviation 
of the perturbative calculation and of the lattice simulations set in. In (d) the instanton cross- 
section in two bins of R/ (p)as the function of the instanton size p is shown. The shaded area 
indicate the experimental limits which is based on the MEPS background description. The 
figures a and b are adapted from Ref. [|7p. 
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